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Abstract 

We present new exact solutions (in 3 + 1 and 2 + 1 dimensions) of relativistic 
wave equations (Klein-Gordon and Dirac) in external electromagnetic fields 
of special form. These fields are combinations of Aharonov-Bohm solenoid 
field and some additional electric and magnetic fields. In particular, as such 
additional fields, we consider longitudinal electric and magnetic fields, some 
crossed fields, and some special non-uniform fields. The solutions obtained can 
be useful to study Aharonov-Bohm effect in the corresponding electromagnetic 
fields. 
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I. INTRODUCTION 



Aharonov-Bohm (AB) effect |T| plays an important role in quantum theory refining the 
status of electromagnetic potentials in this theory. First this effect was discussed in relation 
to a study of interaction between a non-relativistic charged particle and an infinitely long and 
infinitesimally thin magnetic solenoid field (further AB field) (a similar effect was discussed 
earlier by Ehrenberg and Siday 0). It was discovered that particle wave functions vanish 
at the solenoid line. In spite of the fact that the magnetic field vanishes out of the solenoid, 
the phase shift in the wave functions is proportional to the corresponding magnetic flux [[3] . 
A non-trivial particle scattering by the solenoid is interpreted as a possibility for quantum 
particles to "feel" potentials of the corresponding electromagnetic field. Indeed, potentials 
of AB field do not vanish out of the solenoid. AB scattering for spinning particles was 
considered in |4[] and || using exact solutions of Dirac equation in the AB field. A number 
of theoretical works and convinced experiments were done to clarify AB effect and to prove 
its existence (see, for example, 0-0). 

A progress in study of AB effect may be related to revealing new situations, where 
the effect takes place. For example, one can consider more complicated configurations of 
electromagnetic fields, different regimes of particle motions, different dimensions, and so 
on. To study these new possibilities one has to have exact solutions of the corresponding 
quantum equations in these configurations of electromagnetic fields. In this relation, we 
ought to mention exact solutions of the Schrodinger equation in a superposition of AB field 
and a uniform magnetic field fllQf . The latter solutions where analyzed in ||rTHl3|| from AB 



effect point of view. The corresponding coherent states were constructed in |12| . Klein- 
Gordon and Dirac equations for particles moving in a superposition of AB field, Coulomb 
field, and magnetic monopole field were found and analyzed in [p| . |T5|1 . 

In this article we present new exact solutions (in 3+1 and 2+1 dimensions) of relativistic 
wave equations (Klein-Gordon and Dirac) in external electromagnetic fields of special form. 
This fields are combinations of AB field and different types of electric and magnetic fields. In 
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Sect. II we consider AB field combined with longitudinal electromagnetic fields. In Sect. Ill 
superpositions of AB, longitudinal, and crossed fields are studied. In Sect. IV,V we present 
solutions in AB field combined with some non-uniform fields. Here we also discuss some 
relevant solutions in 2 + 1 dimensional QED. Special functions and their properties, which 
are used in the article, are present in the Appendix. 

Most of works, in which AB effect was studied, are based on the use of exact solutions of 
Schrodinger equation in AB field [l]. Consider the latter a field in 3 + 1 dimensions. If the 
magnetic solenoid is placed along the axis z = x 3 , then AB field can be given by potentials 
(we denote these potentials as A^(x), x = (x M , /i = 0, 1, 2, 3)) of the form 

4 0) = 4 0) = ~^ Af = Af = 0, r 2 = (xf + (x 2 ) 2 . (1.1) 

AB mag netic field has the form = (0, 0, H<®), where is singular at r = 0, 

#(°) = ^(x^dix 2 ) . (1.2) 

AB field creates a finite magnetic flux $ along the axis z. It is convenient to define a quantity 
fi, which characterizes the magnetic flux $ and is related to the latter as follows 

$ = (Z + /i )$ , $ = 2-Kchj |e| , < fi < 1, (1.3) 

where Iq is integer, and e = — |e| is the charge of the electron. In what follows, we call 
\i the mantissa of the magnetic flux $. By definition fi is a positive fractional part of the 
magnetic flux if the latter is measured in units of quanta $o- Cylindrical coordinates r, (p 
(x 1 = r cos (f, x 2 = r sin ip) are preferable for AB field consideration. In these coordinates 

= ii±^sin^, = _^ cos ^ . (1.4) 

ch r ch r 

In the present article, we are going to consider particle motion in electromagnetic fields 
that are a combination of AB field and some additional fields with potentials A^ , 

A, = Af + A^. (1.5) 
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Electromagnetic potentials enter in relativistic wave equations only via the operators of mo- 
menta = ihdfj_ — . Doing the transformation ^>{x) = e~ dolfi ^(x) of wave functions, we 
can eliminate 1$ dependence of AB potentials in equations for ty(x). Indeed, such equations 
already contain momentum operators of the form 

e u ovp ^ e -u ov = m ^ _ £^(o) + A (i)^ 

44 0) = -sin^, = -^cos^, 4 0) = if = 0. (1.6) 

ch r ch r 

Thus, all the matrix elements of any axial-symmetric operators depend on the mantissa of 
the magnetic flux only 

II. AHARONOV-BOHM FIELD COMBINED WITH LONGITUDINAL 

ELECTROMAGNETIC FIELDS. 

Here we consider particle motion in a superposition of AB field and of some longitudinal 
electromagnetic fields. We call electric E and magnetic H fields longitudinal ones whenever 
they are parallel and are directed along the AB solenoid (along the axis z). 

E = En, H=Hn, n 2 = 1, n= (0,0, 1) . (2.1) 

It follows from Maxwell equations that in this case the functions E and H must obey the 
conditions 

E = E (x°, x 3 ) = d A 3 - d 3 A , A = A (x°, x 3 ) , A 3 = A 3 (V\ x 3 ) ; 

H = H(x\ x 1 ) = 3 2 A 1 - d x A 2 , A 1 = A l (x\ x 2 ) , A 2 = A 2 (x\x 2 ) , (2.2) 

where Aq (x°,x 3 ) , A\ (x l ,x 2 ) , A 2 (x x ,x 2 ) , A 3 (x°,x 3 ) are arbitrary functions of the indi- 
cated arguments. Exact solutions of the relativistic wave equations in such fields (in the 
absence of AB field) were studied in |T6|JT7|| . As we show below, whenever AB field is 
present, then exact solutions of the relativistic wave equations can be found only in the 
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axially symmetric case with the magnetic field having the form H = H (r) . Thus, po- 
tentials of additional fields, which are considered in the present Section, have the form: 
A$p = ASP (x°, x 3 ) , = A$p (x°, x 3 ) arbitrary and 

,(D chAir) . m chAir) chA'tr) (n „. 

= T-r— -smi/J, A 2 = - — ^^cosy2, if (r) = — r — — , (2.3) 
\e\ r \e\ r \e\ r 

where A (r) is an arbitrary function of r. 



A. Classical description of radial motion. 

To interpret quantum numbers of wave functions, it is often useful to have classical 
picture of the problem. That is why we present here a classical analysis of the particle 
motion in fields under consideration. 

Consider classical trajectories that do not intersect the axis z , thus they do not "feel" 
the existence of AB field. For such trajectories, the quantity P r 2 is an integral of motion 
(c 2 P r 2 is said to be radial energy), 

Pi = Pl + p| = n 2 k\ , mlc 2 + P r 2 = P 2 - P 2 , (2.4) 

where P M is the classical kinetic momentum (a classical analog of the operators P M ) and mo 
is the rest mass. L z is an integral of motion as well (L is the angular momentum), 

L z = l z -h(l + n) =h(l- l ) , L z = x x P 2 - x 2 P 1 - HA (r) =h(l + fi) . (2.5) 

Here I is arbitrary (/ will be integer in quantum theory). 

As will be seen below, exact solutions of relativistic wave equations can be found when- 
ever the functions A (r)in (|2.3|) have the form 

1. A(r) =0, (2.6) 

2 

2. A(r) = T-, 1 >0, (2.7) 

3. A(r) = jr, 7 > . (2.8) 
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The first case corresponds to the absence of an additional electromagnetic field, the second 
one corresponds to the additional constant uniform magnetic field H along the solenoid 
(7 = ^§pj , and the third one corresponds to the additional constant magnetic field H (r) = 
b/r, ( 7 = £A ) . Consider classical motion in these cases. 

1. For A (r) = 0, the momenta Pi and P2 are integrals of motion. Then the radial motion 
(the motion in x^x 2 plane) is parametrized by the proper time r and can be presented as 

pi p2 

x 1 = T + xL, x 2 = T + xL , (2.9) 

where xL^x 2 ^ are integration constants. In this case 

L z = h (I + /i) = xf 0) P 2 - x^P 1 . (2.10) 

Consider the quantity 

AP = P r - 1 (x; 0) P 2 -x 2 0) P 1 ) = ^. (2.11) 



One can show that |AP| characterizes a minimal distance between the trajectory ( |2.9| ) and 
the axis x 3 . All the classical trajectories are divided in two groups according to the sign of 
I + fi . Trajectories with I + fi > can be called right ones and those with I + // < can be 
called left ones. The reason is the following: Looking from the positive z— direction, one can 
see that a minimal angle rotation from the vector r = (x 1 , x 2 , x 3 ) to the particle momentum 
is anticlockwise for the right trajectories and clockwise for left ones. 
2. For A (r) = / jr 2 /2, the radial motion has the form 

x 1 = R cos k + X( ) , x 2 = R sin k + x 2 ) , 

7 TTLqC 

k = UqT + y9 , uj Q = — , m = — — . (2.12) 
m n 

Here R, (p ,x} -,, x? j are integration constants. The trajectories Q2.12j ) are circles of radius 
R with centers having coordinates X( ),x 2 ) , 

(x 1 -xf 0) ) 2 + (x 2 -x 2 0) ) 2 = P 2 . (2.13) 
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One can easily find 

P r = K-yR, I + /I = | (R 2 - Rl) , (xf 0) ) 2 + (x 2 {0) ) 2 = R { 
l + fx < 1 — = -i£- = fi. (2.14) 

We can see that classical trajectories with / > —/i embrace the solenoid, and ones with 
I < —\l do not. In quantum theory these conditions are I > and / < respectively. 
The quantity AR characterizes a minimal distance between the trajectory ( |2.12j ) and the 
solenoid, 

AR=\R-Ro\= j|±^L . (2.15) 
7 (R + R ) 



3. Consider finally A (r) = jr . Here the radial motion depends essentially on values of 
constants a,e, 



P fci 

a = = — > 0, e = sign(Z + a) . (2.16) 

h'j 7 

For e = l, the classical motion is possible only if a > 1. For e = — 1, the classical motion 

is possible if a > 0. Whenever a > 1 we get unbounded motion for r. For < a < 1, 

e = — l,this motion is bounded. Below we present the radial motion in s parametrization 

\l + u\ (a cosh s + e) \l + u\ (a sinh s + es) 

a > 1 : r = , : , t - 



7(« 2 -l) ' " 7 2 (« 2 -l) 3/2 



g f d — £ g 

ip — loq = -. + 2e arctan W tanh - 

Va -1 \ V a + £ 2 / 

a = 1, £ = -1 : 27r = \l + /i| (s 2 + l) , 2 7 2 r = |Z + /i| m + s^J 

99 — ip = s — 2 arctan s ; 

1/ + u\ (1 — acos s) \l + u\ m (s — a sin s) 

a < 1, e = -1 : r = J ^ - -, r~ 



7(l-« 2 ) ' ~ 7 2 (l-a 2 ) 3/2 



/ 1 \ / a sins \ 

V ~ ¥0 = si R -1-2 arctan 7 = . (2.17) 

\V1 - fl / \ 1 + v 1 — a —a cos sj 

In all the cases under consideration, the minimal distance between a trajectory and the 
solenoid is defined by the expression 



Thus, the quantity / has a clear classical interpretation. 



B. Klein-Gordon equation in longitudinal fields 

Here we consider solutions of Klein-Gordon equation 

P M P M - m 2 c 2 ) * (x) = (2.19) 

in the superposition of the external fields ( |1.1| ) and fl2.3|). In this case, the operators ( |2.4|) 
are integrals of motion. Whenever an additional field is axial-symmetric one ( |2.3| ), then the 
operator ( |2.5| ) is an integral of motion as well. Thus, we subject solutions of the equation 
( |2.19j ) to the following additional conditions 

n -2 f p 2 + p2\ ^ ^ = k 2y ^ ^ n -2 f p 2 _ p2\ ^ ^ = ^2 + ^ ^ _ ^ ^ 



Then such solutions can be presented in the form ^ (x) = ip (x\ x 2 ) <3? (x°, x 3 ) ,where the 
functions ip and $ obey the equations 

h- 2 (p 2 + P 2 )^{x\x 2 ) =k\i,(x\x 2 ) , (2.21) 
h- 2 (P 2 - P 2 ) $ x 3 ) = (m 2 + kl)<$> (x°, x 3 ) . (2.22) 

AB field does not enter in the equation ( |2.22| ). This equation can be solved for a large 
class of electromagnetic fields. All the corresponding solutions of the equation ( |2.22j ) are 
described in detail in [ IB| , pT |, that is why we do not present them here. 

Let us turn to the equation ( |2.21| ). As was already mentioned above, exact solutions of 
this equation can be found only in the superposition of AB field and the fields fl2.6|)-(f2~8|). 
In all these cases, the operator L z is an integral of motion, thus we can search for solutions 
which are eigenvectors for the latter operator. In cylindrical coordinates, we get 

L z ip (r, <p) = -ihd^ (r, <p) = h(l- la) * (r, <p) , V (r, <p) = ^M^M ^ ( r ) . (2 .23) 

V 27T 



Whenever additional fields have the structure (|2.3|) , the radial function ip (r) obeys the 
equation 

/ / + /i + A(r)\ 2 Id d 2 



Rip (r) = kfijj (r) , R 



r dr dr 2 



(2.24) 



C. Solutions of radial equation in the absence of additional fields 



Consider the equation (|2.24|) for A (r) = 



d 2 ip 1 dip 
dr 2 r dr 
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if) = 0. 



(2.25) 



A general solution of the equation can be written via Bessel functions J u (x) , ( |]I8"||, 8.402) 



i> ( r ) = ( r ) = (^i r ) + c 2 ^- y , v = \l + n\ . 



(2.26) 



Solutions, which are bounded for all r > 0, must have Ci = 0. In this case ipki,i form an 
orthogonal and complete set of functions (ki, k[ > 0), 

poo poo () (r v ,s \ 

/ J v (k[r) J u (hr) rdr = jferM (&i - tf) , / J u {hr') J u {hr) hdh = -i ^. (2.27) 

Jo jo r 

Suppose Z = 0, —1, fi 7^ (note that in this case the corresponding classical trajectories 
pass maximally close to the solenoid). Then we have a special case. Here there are solutions 
of the form ( |2.26p with c\ = 0, c% ^ (of course they are unbounded), which obey the 
relations ( |2.27|) . Moreover, there are solutions (|2.26|) with both c\ ^ and Ci ^ 0, which 
are unbounded and not orthogonal, 



J v (k[r) J_ u {hr) rdr= , 
(I 7i (k± — kf) \k\ 



2 sin z/7r /^V , 
— — t^ItH > *1>*1, H<1- 



(2.28) 



However, these solutions are quadratically integrable (due to (|2.27|) ) and form a complete 
set of functions (in fact, an overcomplete set). Besides, the equation ( |2.25| ) has solutions of 
the form 



■0 (r) = K v {qr) , k\ 



(2.29) 



where K v (r) are Macdonald functions ( ||18|| , 8.407). These functions have a finite norm 
(5 = arg q) 

roo 7j- c[n 2uS 

/ K* (qr) K v (qr) rdr = - 2 , \v\ < 1 , -vr/2 < 5 < tt/2 (Reg > 0). (2.30) 

Jo 2\q\ sin z/7r sin 26 

But they are not orthogonal with respect to q, 

. K v (q'r) K v (qr) rdr = - 2 {q \ ) , \u\ < 1 . (2.31) 
o 2{q 2 — g /z J sin utt 



In particular, for v = 0, we get 

/ ^(gr)K (gr)rrfr= / K (q'r) K (qr) rdr = . (2.32) 

Jo \a\ sin 2o Jo Q — (7 



The above peculiarities are related to the loss of hermicity of the operator R for I = 0, — 1, 
/x 0. Similar problem was discussed in |19| . 



D. Uniform magnetic field 

Here we consider a superposition of AB field ( |1.1| ) and a uniform magnetic field ( |2.7| ). It 
is useful to introduce dimensionless operators a&, ajj", k = 1, 2 by the relations, 

fty^ = -iP a - P 2 , fry^a^ = zP x - P 2 , 7 = ^ , 

ft^as = -zPi + P 2 + h^(x l + ix 2 ), h^aj = iP 1 + P 2 + h^(x l - ix 2 ) . (2.33) 



Considering coordinates and momenta in these relations as classical quantities, we can get 
a representation for the classical motion ( 2.12j ) in terms of a\ and a 2 



a x = ]f^Re~ lK , a 2 = ^l{x\ + ixl) , {x\ + ix 2 Q ) = R e lS . (2.34) 
The following operator relations take place 

P 2 = P 2 + P 2 2 = fi 2 ^(a^ai + aiaf) , 2L Z = h(afai + a\a^ — a^a 2 — a 2 a 2 f ). (2.35) 
We introduce also a dimensionless coordinate p instead of r, 
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'yv 2 1 

p = , dx l dx 2 = —dpdip . (2.36) 

2 7 



On the classical trajectories fl2.12| ) p evolves as 

2p = 7 \R 2 + R 2 + 2RR cos (k - 5)] . (2.37) 

Being written in terms of the variables p, ip, the operators a&, al take the form 

ai = y/pe-'v [(l + p + p- id v )/2p + d p ] , 
4 = 4~P^ Wo + P + P~ id v )/2p - d p ] , 
a 2 = ~y/pe iLp [(k + fi- p- id v )/2p - d„] , 

at = ~VP e ' i<P i( l o + P-P- id 9 )/2p + d p ] . (2.38) 

Using the commutation relations for the momentum operators 

eh 

PpP v - P v Pp = -i—F^, Fp, v = dpA v - d u Ap , (2.39) 



and definition of the magnetic field (|1.2|) , we arrive to the following commutation relations 
for the operators a,k, 

[ai,at] = l + f, [a 2 , ajj"] = 1 - /, [a x , a 2 ] = -/, [ai,a^]=0, 

/ = (^/H)5(x 1 )5(x 2 ) = 2§-5(p) =2(l + p) 5(p) . (2.40) 

These commutation relations contain a singular dimensionless function /. Whenever AB field 
is absent (/ = 0), then the operators ak,a^ form two mutual commuting sets of creation 
and annihilation operators. It is not true in the presence of AB field. However, as it will 
be seen further, these operators behave as creation and annihilation ones when acting on 
functions that tend to zero (sufficiently rapidly) as p — > 0. Being written in the coordinates 
p, ip, the operators ( |2.35|) have the form 

P 2 = 2~fh 2 Q , L z = -ihd v , Q = ik + fl+ / p ~ %dip)2 -d P - P d 2 p . (2.41) 

In the case under consideration, the radial equation Q2.24 ) reads 

11 



W(P) = (n + \)^(p) ,Q= il + ^ P)2 -d p -pd 2 p , kf = 2 1 (n+^J . (2.42) 

Bounded and quadratically integrable solutions of this equation are expressed via the La- 
guerre functions ( |6.1| ) (see Appendix). There are two types of solutions the latter equation, 
we denote them as ijj^'(p), j = 1,2 (two types of states). The first one j = 1 corresponds 
to / > (classical trajectories embrace the solenoid) 

V (1) (P) = 0<l<n, n = 0,1,2,..., n = n + /i . (2.43) 

The second type of solutions with j = 2 corresponds to / < (classical trajectories do not 
embrace the solenoid) 

iP {2 \ P ) = In-i-„,n(p) , 1 < °> n = n- (2.44) 
In these two cases radial momentum spectra are different, 
(k^Y = 2j (n + ^ + , < I < n, 

(k?) 2 = 2 7 (n+^,l<0, n = 0,l,2,.... (2.45) 

The spectrum for j = 2 (which is a part of the total spectrum) corresponds exactly to the 
spectrum of a spinless particle in a uniform magnetic field (without AB field). The spectrum 
for j = 1 is shifted by 27/x with respect to the one for j = 2. It is important to note that 
the presence of AB field lifts partially the degeneracy of the total spectrum in the quantum 
number I. 

It is convenient to define effective quantum numbers / and n by the relations 

in + u , 7 = 1, I = I + u, I < n 
(2.46) 
n, j = 2, n = 0,1,2,.... 

Using these numbers, we introduce the functions 

/(l)/ \ / ^n-i ex P W ~ l M r 1 \ 

*Pn l(P, V) = 7== 4,n-KP) , 

V 271 



€}M={ _ lT ^m^m h _ lM . (2 .4 7 ) 

V 2n 
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According to ( |6.4| ), these functions can be expressed via the Laguerre polynomials. Thus, 
the orthonormality relation can be proved 



oo p2tt . 

dp / dip ipzjr (p, (p)ip%j (p, <p) = 6 Ll >6 n , r 
o Jo 



(2.48) 



The set of the Laguerre functions 

I a +n,n{x), n = 0, 1, 2... , a > -1 
is complete in the space of quadratically integrable functions of x > 0, 

oo 

Ia+n,n{x)Ia+n,n{y) = 5 (x - y) . 

n=0 



(2.49) 



(2.50) 



Then the set ip^ ] (p, (p) is complete in the space of quadratically integrable functions of p, (p, 
(p > 0, < <p < 2tt). 

Using the relations (|6.6|)-( |6TTTD , one can get the action of the operators ( [2.3S| ) on the 
functions xp^](p,(p) 7 



OiV>nJ(/w) = vWn-i,/-i(p,^), a f^n}(p^) = Vn + l^i+i, m (p, (p), 
a2ipn}{p,(f) = J™ ~ Wnj+iiP, <P), at^n}(P,f) = "Jn - I + lipnj^ip, (p). 



(2.51) 



These formulas show that the functions ip^ } ma Y be created by an action of the operators a] 
on t/>qo , and the functions i/j^j may be created by an action of the operators a~l on 



Namely, 



(i) 



T(l + p) 



\ r(l + n)r(l + n - o 



(4r l (atMl 



T(2-p) 



\ r(i + n)r(l + n - o 



[2.52) 
[2.53) 



It is natural to interpret as a vacuum state for the states \ ■> an d to interpret ^o 2 -i 

(2) 

as a vacuum state for the states ip n \ . Thus, for p ^ 0, we have two vacuum states in the 
problem. For p = 0, the situations changes. By virtue of ( |6.18 ) 



'n,n— I 



-I) 1 I n ^ n -> ^} 



(-i)VS. m = o 



(2.54) 
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and for any I < n, the function t/>qq is connected to V'cj 2 -! as 



„+./.(!) _ l/; ( 2 ) „ 7 /,(2) _ „/.(i) 

a 2 ^0,0 — ¥>0,-l> a 2Wo -1 — ^0,0 • 



(2.55) 



Thus, we have the only one vacuum in the problem, one energy spectrum ( |2.45|) , and all the 
wave functions are created from the vacuum V'oo- 

One ought to stress that all the states obey the property ip„] (p — 0, (p) = , which 
means that the scalar particle has zero probability to be found in the solenoid area. In fact, 
the existence of this property allows us to speak about AB effect. 

The definitions ( 2.51) can formally be considered for any values of indices n, I. In par- 
ticular, we can consider the following relations 



■l) n+1 (l+n) 



ai#o = V n + 



<" + "H arrpVVln) exp| " !(1 + h)v - l ]/, ~* L "- 1M ' 



a>ti>n,o = v / n + W n ,-i 

sn+l I 



n 



r(l + n) r , . p, i-M . , N 



(2) 

1 - p + ra/'n.o 



■lf(l-p + n) 



However, the functions ^i-i; V'no do no ^ P re sent any physical solutions of the problem, 
they are not in the set ( [2.47 ). Thus, in the general case p 7^ 0, the action of the operators 
afc, cik on wave functions may lead them out of a class of physical solutions. The functions 
( p.56| ) are singular at r = (for p 7^ 0), however they still remain quadratically integrable. 

Thus, we see that I = 0, —1 is a special case. Here there appear unbounded (but 
quadratically integrable) solution V'i-iJ ^nl- Whenever p — > 0, these states either coin- 
cide with the corresponding states in the pure magnetic field or disappear. The states 



^ 27 rr(2-p + n) eXp[ -^-2 ]P " L - ^ 



(2.56) 
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^io; V'nO) V'n-i) Vn 2 -i are n °t mutually orthogonal in spite of the fact that they belong to 
different eigenvalues of the operator P 2 . 

The equation ( |2.42|) has additional solutions in the case I = 0,-1. According to (p. 20 ) 
they have the form 



4>(p) = *Pxa(p), a= { ' , 2n = 2X + l + n-l, (2.57) 




where the functions i[)\ a (p) are defined by Eqs. (|6.28|) , (|6.29|) . Solutions ( 2.57|) with any 



different complex A are orthogonal and have finite norms according to the properties 



( p. fC ). These solutions are singular at r = 0. For I = such solutions exist even in the pure 
magnetic field. Their existence is related to the loss of hermicity of the operator P 2 . 

E. Nonuniform magnetic field 



Here we consider the radial equation (3.2.8) for A (r) = jr, 

|t + 1< T ~ ( + * = 0, P; = (n,af , kl = (;af . (2.58) 
ar z r ar \ r 



The constant a is defined in Q2.16| ). 



For a 7^ 1 bounded solutions of this equation are expressed via the the Laguerre functions 
.ll) as 



i/> (r) = ipnj (r) = I a +n,n(x) , x = 2\J \ - a 2 ^r , 

2(1 + u) , 

a = 2\l + fi\ , l + a + 2n = )==■ (2.59) 

V 1 — a 2 

For a > 1, there are solutions for any I (in complete accordance with the classical theory). 
In this case the Laguerre functions have imaginary arguments and complex indices. 

For a < 1, bounded solutions of the form (|2.59| ) exist only for / < (also in accordance 
with the classical theory). Besides, in such a case 

a 2 2aryr 

a 2 = 1 ? , n = 0,1,2, x= [ . 2.60 

(l + a + 2nf l + a + 2n y 1 
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Thus, n must be integer and a is quantized. Here the functions (|2.59|) can be expressed via 
the Laguerre polynomials by means of ( |6.4j ), ( |6.19| ). 



For a = 1, bounded solutions of the form ( |2.59| ) exist only for I < 0, these functions can 
be expressed via the Bessel functions, 

V> (r) = J a (2-y/ayf) , a = 1, I < . (2.61) 



One can see with the help of Q6.22D that the solutions ( 2.61|) follow from ( 2.59 ) as a — > 1. 



All the bounded solutions vanish at r = 0. 

It is interesting to note that there exist unbounded (but quadratically integrable solu- 
tions) of the equation ( |2.58| ) for I = 0,-1 . For any complex a, the latter solutions are 
defined as 

i= (r) = ipx, a (2v / 1^70 , a ^ 1, A = , Re v / T^ > , 

v 7 v 1 — cr 

^ (r) = K Q {2yfcyP) , a = 1 . (2.62) 

Here (x) are Macdonald functions. The existence of such solutions is related to the loss 
of hermicity of the operator P r 2 . 



F. Solutions of Klein-Gordon equation that are not related to radial momentum 

conservation 

As was demonstrated above, selecting the radial momentum ( |2.21|) as an integral of 
motion, we can separate variables and then consider two independent problems: a two- 
dimensional motion of the charge in the mag netic field (|2~2l| ) (the latter field includes AB 



field), and a two-dimensional motion of the charge in an electric field, the latter problem 
does not depend on AB field. However, there is a wide class of exact solutions, which are not 
eigenvectors for the radial momentum operator. They correspond to a superposition of AB 
field and longitudinal running electric fields (potentials of such fields depend on u° = x° — x 3 
only). Thus, here we will use light cone variables u°,u 3 , 



16 



u° = x° - x 3 , u 3 =x° + x 3 



(2.63) 



Then the above mentioned longitudinal running electric fields have the following potentials 
and strengths 



1 



B(u u ), E = B' (u° 



A {1) - A {1) 

where B (u°) is an arbitrary function of u°. Consider operators P ,P 3 ,p 3 





(2.64) 



2P =Po- Pz, 2P 3 = Po + Ps, Po = ih 



- -fi 9 P -a. ^ 



e| B (u°) 

ch 



[2.65) 



The operator p 3 commutes with the one L z and both are integrals of motion. Thus, we can 
demand for solutions of Eq. ( [2.19| ) to be eigenvectors for these operators, 



hX 

p 3 y (x) = — V (x) , L z m (x) =h(l-lo)* (x) . 



(2.66) 



Such solutions have the form 



iff (x) = X + g 1/2 $ (r, t) exp i 

du° 



(I -l Q )(p- mH (u°) - 



Xu 3 



t 



X + g(u°)> 

where the function $ (r, t) obeys the equation 



(2.67) 



Ri$(r,t) = 0, R^idt + d- 



Or 



l + A(r) 



1 = 1 + fj, . 



(2.68) 



We recall that A (r) was defined in (|2[ 

Consider first the case A (r) = 0. Here we find a propagation function for the equation 
( |2.68 ) in the form (J u (x) are the Bessel functions) 



^.2 j ^'2 



R 1 a=q Go (r, r', t) = , Urn G (r, r' , t) = -8 (r - r') . 



1 7T 



(2.69) 
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The case A(r) = p = / yr 2 /2 can be considered in the same manner. Here the propagation 
function has the form 



P + P' 



G(p,f/,t) = (^) e*>, Q 

zsmr I I V smr J 

R 1 G(p,p',t) = 0, KmG(p,p',t) = 6(p-p') , r = jt (u°) . 



It. 



smr 



(2.70) 



The functions Go { r , r> , t) and G (p, p', t) solve the Cauchy problem. For example, 

$(p,t)= / G(p,p',t)$(p')dp', (2.71) 



where $ (p) is an arbitrary functions (an initial date for $ (p,t)). 

For the field ([2.8|), the corresponding propagation function is quite complicated pO . 



G. Exact solutions of Dirac equation 



Here we are going to study Dirac equation 



(7"P„-moc)tf (a;) =0 



(2.72) 



in the superposition of AB field and the field ( |2 . 1| ) . We use a standard representation (see 
for example []I7|| ) for 7— matrices. In the case under consideration, we look for solutions with 
a definite radial momentum. The corresponding bispinors ^ (x) can be written in a block 
form 



# (x) = Q 



f ip 1 (x 1 , x 2 ) [m + F — ikia 2 ] ^ 



v$ (x°, x 3 ) , F = h- 1 (P + Pa) , (2.73) 



. (x 1 , x 2 ) [(m - F) o- 3 - iA;i<ri] , 

where v is an arbitrary spinor; (k = 1,2,3) are Pauli matrices; the function $(x°,x 3 ) 
obeys the equation 

\e\E~ 



ch 



<§ (a; , a; 3 ) = (m 2 + k\)§ (. 



x°,a; 3 



^2.74) 



where E = E (x°, x 3 ) is electric field strength ( |2.2| ), and functions ip2 obey the following 
equations 



(Pi + iP 2 ) V>i (x\ x 2 ) = hfa^ 2 (x 1 , x 2 ) , (P - iP 2 ) i) 2 (x\x 2 ) = hfa^ (x\ x 2 ) . (2.75) 



The presence of the arbitrary spinor v in the solutions (|2.73|) indicates that Eq. (|2.72|) does 



not fix the spin orientation. This orientation can be fixed by a choice of a spin operator 
1^ , 17] . One has to stress that the equation (|2.74j) does not contain AB field. All possible 



exact solutions of this equation were presented in ||16| , |17[ , thus here we do not repeat these 
results. 

The fields ( |2.3| ) are axially symmetric, thus J z is an integral of motion in such a case 
(J is total angular momentum operator). Let us consider solutions that are eigenvectors for 
this operator, 

J,* = h (l - l - ^) * , J z = L z + -E 3 , 1 = 0, ±1, ±2, ... . (2.76) 
(S =diag (a, a) ).We obey the equations ( |2.76| ) choosing 

V / -\ / /IT \ / / / tt 



where the functions ( r ) satisfy a set of first order differential equations 

( l -±^± + ±) i, 2 (r) = fa^ (r) , P~ 1+A(r) - ^ (r) = jfc^ (r) . (2.78) 
\ r ar J \ r ar J 

Consider solutions of the latter equations for the fields (|2.6|) -(|2~8|). 

For A (r) = 0, we deal with the pure AB field. For I ^ 0, all bounded solutions of Eq. 

( [2.781 ) have the form 

V>i (r) = Ji^i (Atir) , ^ 2 (r) = eJ|,-j (for) , £ = sign / , (2.79) 

where J M (x) are the Bessel functions. These solutions vanish at r = 0. For 1 = 0, \x ^ 0, 
the system of equations ( |2.78| ) has no bounded solutions. In such general solution 

of this system has the form 

ipi (r) = ci J M _i (for) + c 2 Ji_ M (for) , ip 2 (r) = ci (far) - c 2 J-^ (far) . (2.80) 

where ci, c 2 are arbitrary constants. In spite of the fact that these solutions are unbounded 
they still are quadratically integrable (as in the scalar case). Moreover, for any complex fa 
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(Re ki > 0) there exist unbounded solutions with a finite norm, they are expressed via the 
Macdonald functions, 

V>i (r) = K x ^ (hr) , i) 2 (r) = -K^ (hr) , < p < 1 . (2.81) 

Similar to the scalar case, we can conclude that the operator 7 M P M is not self-conjugate 
anymore for I = 0, p ^ 0. In contrast to the scalar case, there are no quadratically integrable 
unbounded solutions for I — — 1, as well as for I = 0, p = 0. 

Consider now the case of the uniform magnetic field ( |2.7| ). Using the operators ( f2.38| ), 
we can write the equations ( p . T5|) as 

2 

ai4>2 (p, <p) = -iyMipi (p, (p) , aft/i! (p, <p) = iVni/j 2 (p,<p) , h = ^2>yn, p = — . 

Their solutions have the form (see Q2.51Q , ( p.47|) ) 

^ (p, ip) = ^i,,.! (p, f) , ^2 (p, if) = (p, f) ■ (2.82) 

As in the scalar case, there are two types of states (with j = 1,2). These states are bounded 
at I 7^ 0; they vanish at r = 0. The states (|2.82|) are unbounded at I = but they still are 
quadratically integrable. Besides, there are unbounded solutions with finite norms for any 
complex n. Such solutions are expressed via the functions ijj\ tCX (x) (the latter are defined by 
( PI ), (PI) as 

^ ( r ) = n 3 / 4 VA-i,i-M (p) , ^2 (r) = n 1 / 4 ^,,, (p) , 2n = 2A + p - 1 . (2.83) 

Thus, we see that the operator 7 M P M is not self-conjugate for I = 0, p ^ as well. All the 
above singular solutions vanish or become nonsingular as p — > 0. 

One ought also remark that ?/>i (which correspond to j = 2) from ( |2.82[ ) vanish at n — 0. 



Thus, the complete wave function Q2.73j ) is an eigenvector for the operator E 3 , 

S 3 ^ = -^. (2.84) 

That means that in such states the electron spin has the only one orientation, namely, 
opposite to the magnetic field. 
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Consider finally the case of nonuniform magnetic field fl2.8|) . For I ^ 0, a ^ 1, the 
corresponding bounded solutions (they also vanish at r = 0) have the form 

ipi (r) = 4_i >n+ i_2/ (x) , V>2 (r) = -I n , n - 2 i(x) , I > ; 



^i (r) 



z n+l-2« 



,n-l (x) , ^2 (r) = -I n -2l,n (x) , I < 



i _ 1-21 - 

x = 2V1 - a V , 1 - 2/ + 2n = , 1 = 1 + {i. 



(2.85) 



where J nj . m (x) are the Laguerre functions (H7TT) , and we use the notation (ETTBI). Whenever 
a 2 > 1, any / ^ are admissible in complete agreement with classical theory. Whenever 



< 1, the only I < are admissible. In such a case n is integer and the functions Q2.85 ) are 



expressed via the Laguerre polynomials according to (|6.19| ). At the same time, the following 
quantization takes place 



1 + 2 



1 + 2 



2 , n= 0,1,2, 



(2.86) 



+ 2nj 

For a = 1, I 0, the only bounded states can be found for I < 0. They are expressed via 
the Bessel functions, 



(r) = J 2|r|+2 ^(l + 2|j|) 7 r) , ^2 (r) = - ^ (l + 2 |/|) 7 r) . 



(2.87) 



Solutions ( p.87| ) follow from ( |2.85|) as a — > 1. That fact can be confirmed by the use of the 
limit ( ggg ). 

Z = is a special case. Here there are only unbounded solutions. Some of them are 
quadratically integrable. Whenever a 2 > 1, such solutions have the form 



if>! (r) = Ci/ n+ i_2 M ,n-i (a) + c 2 4-i,n+i-2 A1 (a) , 

^2 (r) = -Ci/ n _ 2Min (x) - C 2 I n ,n-2t, (x) , 

where are arbitrary constants, and for a = 1 these solutions read 

ipi (r) = J 2 -2 fJ , (2) , ^2 (r) = - J-2/j, (z) , < /i < ^ ; 

V>i (r) = iT 2 -2^ (2) , ^2 (r) = -ftT 2M (2) , - < \i < 1; ^ = 2 v /|l-2/i| 7 r 
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(2.88) 



(2.89) 



Quadratically integrable solutions exist for a 2 < 1 as well. For example, for < // < | they 
have the form ( |2.88|) , where C2 = 0. In such a case a 2 is quantized 



1 - 



l-2/i + 2n) 



2 , n= 0,1,2, 



(2.90) 



Moreover, for any complex a 2 (provided Re \J\ — a 2 > 0) there exist unbounded solutions 
with a finite norm. They read 

1-2/i 



V>i (r) = a^A,2(i-ft) (a:) , ^2 (r) = (l + Vl - a 2 ) ip x ,2^ (x) , A 



2s/T^l? 



[2.91) 



All the above solutions obey (|2.84j ) for n — 0. 

Finally we present solutions, which do not have an analog in the Klein-Gordon case dis- 
cussed in Sect. II. F. These solutions are not eigenvectors of the radial momentum operator. 
With this aim in view we present Dirac wave functions in the following form 



* (x) = tt H Or) + ^ {+) (x) , * (±) (x) = P (±) * (x) , 2P {±) = 1 ± (an) 



(2.92) 



where n is a unit vector, a = (ak = 7°7 fc ) , k = 1,2,3, and P(±) are projection 
operators, P( + ) + P(_) = 1 , P 2 ± ) = P(±) , P( + )P(_) = P(_)P( + ) = . Then we can always 
present \P(±) (x) in the following block form 



*(+) W 



v (x) 
v (trn) v (x) y 



*(-) (*) 



u (x) 
\ - (^n) « (a?) y 



(2.93) 



with w (x) , v (x) being arbitrary spinors. Without loss of generality we can always chose 
n = (0, 0, 1) . The Dirac equation ( |2.72|) demands u (x) , v (x) to obey the equations 



2P m = [m c — (crB) 03] v , 2P 3 -u = [m c + (crB) cr 3 ] u . 
2P = Po — -P3 , 2P 3 = P + P 3 , B = (P 1; P 2 , 0) . 



(2.94) 



Suppose we consider external fields, for which the operator P3 ( |2.65| ) is an integral of motion, 
and suppose we are looking for solutions that are eigenvectors of the latter operator. Then 
in accordance with ( |2.66 ) 
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2P 3 = h(\ + g), chg=\e\ (A™ + A™) 



(2.95) 



It follows from fl2J| ) that the spinor v can be restored by the one 



it, 



h (A + g) v = [m c + (crB) 03] u . 



(2.96) 



For external fields under consideration, the operator B commutes with A + g, thus we get a 
closed equation for u, 



2h (A + g) P u = [m c — (crB) <7 3 ] [m c + (crB) cr 3 ] w . 



(2.97) 



Considering eigenvectors for the operator J z (|2.76|) in axial-symmetric external fields Q2.3|) , 
we can write 



' e-^ui (r, t) » 



MX 



exp z 



(/ - Z )<£ 



Am 3 



where the functions (r, t) , ( = ±1 obey the equations 



2 7 A + #(w°) 2 



(2.98) 



R{u c (r, t) = , = iSt + d 2 + — - 



l- l -^f + A{r) 



-c 



A' (r) 



(2.99) 



which can be solved similar to the one Q2.68 ). 



III. SUPERPOSITION OF AHARONOV-BOHM, LONGITUDINAL, AND 

CROSSED FIELDS 



We consider here Klein-Gordon and Dirac equations in some superpositions of AB field, 
longitudinal, and crossed fields. In fact, there are only two types of such fields, which admit 
exact solutions of these equations. 

To define the first type of the fields, we introduce curvilinear coordinates by the 
relations 



u° = x° - x 3 , u 



1 = q (u°) r 2 , u 2 = (p, u 3 = x° + x 3 - u%\ q (u°) = Uu )" + 



(3.1) 
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where a is a constant. In these coordinates, covariant components A\p of electromagnetic 



potentials are given as 



h4" 

ch 



q(u°) [h(u')+au l gi (u )}, A? = , 



f2(u 1 )+u 1 g 2 (u°) , 



\e\A 



(i) 



9i (u°) 



(3.2) 



eft ^ v / J * v y 2 

Here g s (w°) , / s (u 1 ) (s = 1, 2) are arbitrary functions of indicated arguments. The cor- 
responding additional (to AB field) electromagnetic field is given by its components in 
cylindrical reference frame 

\e\E z 



e\E r 


\e\H<p 


ch 


ch 


e\E v _ 


\e\H r 


ch 


ch 



qr 



2g(/{ + agi) + u°g[\ , 
qr[g' 2 -2qu°(g 2 + f^}, 



ch 

\e\H z 

ch 



~9i, 

= 2q{g 2 + f 2 ). 



(3.3) 



Exact solutions in the field ( |3.3| ) were studied in [f21] ,|T7jl. 

In the case under considerations, the operators L z (or J z in the Dirac equation case) and 



P3 ( |2.65| ) are integrals of motion. We are going to study solutions that are eigenvectors 
for such operators. Let us impose the following constraint on the functions g s (u°) , f s (u 1 ) 
(a = 1,2) 



u 1 (gl - ag{ -b) + 2g 2 f 2 - 2g 1 f 1 = 0, b = const . 



(3.4) 



Then we can separate the variables u° and u 1 and present Klein-Gordon wave functions in 
the form 



q (u 1 ) , P = X + g x (u 1 



X f dvP 

T = -u 3 -(l-l )<p + J [m 2 + 2q (2 A* + Ig^ ' ' 



2P 



(3.5) 



The functions ip (u 1 ) obey the equation 



ip" + + R (u 1 ) tj) = , R [u 1 



2fci + A/i aX 2 + b [l + f2 



2u l 



(3.6) 



In the same case, Dirac wave functions have the form 

24 



p 



where v is an arbitrary constant spinor, and 



(l + <7 3 )Vl M +(1-^3)^-1 



(3.7) 



7 m + P-a 3 (aF) ^ 



F = e r gr (2i<9 u i - w°P 



/ 



v (m - P) a 3 - (aF) 

Z_ +/2 . \ TTr ... ... /"g^2 , 

h qrg 2 , W = cos + ta 3 sm 0, = aw . 

r J P 



(3.8) 



The scalar functions ipc ( ul ) (C = ±1) obey a set of independent equations 



4' + -^c + 



(3.9) 



Explicit solutions of the equations ( |3.6| ) and ( |3.9| ) can be written for 
fx (u 1 ^ = au 1 + -^j- , / 2 = 7 = const 



(3.10) 



In such a case the equations (|3.6| ) and (|3.9|) are reduced to the one ( [2.581 ) . Solutions of the 
latter equation we have studied above. 

Let us return to the constraint (|3.4| ). If (3 ^ 0, then gi = 0, (72 = const, and b can be 
found from (|3~4]) to be b = g\ + 27(72 . If (5 = 0, then g>i, g 2 are related by an equation 



(5-2 + af = a (gi + + 



2 <-y 2 

2 / 1 7 

a ho. 

a 



(3.11) 



Thus, one of the constants remains arbitrary. We see that there exist a wide class of fields, 
which admit exact solutions. 

To define the second type of fields, which admit exact solutions, we introduce curvilinear 
coordinates by the relations 



u° = x° - a; 3 , u 1 = —,«" = <f, u° = x" + x° . 



(3.12) 



Covariant components of electromagnetic potentials in the coordinates ( |3.12| ) are given by 
the equations 



|e|4" 

ch 



h (u 1 ) 



A? = 0, 



|e|4 1} 

ch 



A W 

j ^3 



(3.13) 
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where f s (u 1 ) (s = 1,2) are arbitrary functions of u 1 . The corresponding electromagnetic 
field is given by its components in the cylindrical reference frame 



E r = H,, 



2chr 



2chr 



2ch 



/£, H z = —f 2 , E z = 0. (3.14) 



e u 



|e| [u u y \e\ (u°) 

In the absence of AB field, exact solutions in such a field were studied in [^2 ,17 



Here integrals of motion are the same as in the previous case. Klein-Gordon wave func- 
tions can be written in the form 



^u 3 — i(l — l )ip + m 2 u° + ki lnw 



(3.15) 



The functions ip (w 1 ) obey the equation 



u v 1 v J 16 4m 1 4 (u l 



(3.16) 



16 4m 1 4(u x ) z 
Dirac wave function can be presented in the form ( |3.7| ), ( |3.8| ) with the following modifications 

l + h 



1 r A 

P = X, q=— ,F = e r — [ 2z9a 



M° ' M° V""* 2 



(3.17) 



Besides, the functions ^ (w 1 ) have to obey the equation ( |3.9|) with R (u 1 ) defined by ( |3.16|) . 
Solutions of the latter equations are available for f s (u 1 ) in the form ( |3.10| ). Thus, we have 
again returned to the equation (|2.58|). 



IV. SUPERPOSITION OF AHARONOV-BOHM FIELD AND SOME 

NON-UNIFORM FIELDS 



Consider now additional fields, which are given by potentials of the form 



4 1} = n/i(r), A 



ch 



(i) 



ch A (r) 



sin <p , A 2 



(i) ch A (r) 



cos<^, A^ = ^f 2 (r) . (4.1) 



ch 



\e\ r \e\ r \e\ 

Here /i (r) , f 2 (r) , A (r) are arbitrary functions of r. The corresponding electromagnetic 
field components in cylindrical reference frame have the form 



ft = -n/'.M,ir, = n/.'M,ff. = n— 

e e e r 



, E v = E z = H r = . 



(4.2) 
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Exact solutions of relativistic wave equations in such fields were studied in p3] , |T7|l . Here we 
present exact solutions of the equations in the superposition of these fields and AB field. 

Stationary solutions of Klein-Gordon equation that are eigenvectors for the operators 
Po,P3, L z can be written as 



$ ( x ) = e - ir ^> (r) , T = k x° + hx A - (l-l )<p. 



(4.3) 



Functions tp (r) obey the equation 



( r ) + Itf ( r ) + R (r) $ ( r ) + o , R (r) = {k + hf - (k 3 + / 2 ) 2 - ^ V " m * ■ ( 44 ) 

The corresponding solutions of Dirac equation have the form \1/ (x) = e~ iT 'Mi/j, where the 
matrix M reads M = diag(e~ lip ,i,e~ lip ,i) , and the bispinor ip=(ipk) (k = 1,2,3,4) obeys 
the equation 



ka + fi- m 7 ° - i (l + A-~^ ai -i\jt + ^ja 2 + (h + f 2 )a V = . 



(4.5) 



In some particular cases the latter equation can be reduced to the one ( |2.78| ) and thus solved 
explicitly. All such cases are described in P3|JT7| . 

V. AHARONOV-BOHM FIELD IN 2+1 QED 



Consider Dirac equation in 2 + 1 dimensions (see for example, [24,p5[) (x = (x M ), \i 



0,1,2, 7° = a 3 , 7 1 = ^ 2 , 7 3 



-zov 



( 7 M P M - m c) * (x) = , * (x) 



#2 (x) 



(5.1) 



/ 



For components of the spinor \1/ (x) we get the following equations 



(P - m) *i + (Pi - iP 2 ) ^2 = 0, (P + m) ^ 2 + (Pi + 1P2) * 1 = 



(5.2) 



These equations can be solved exactly for a superposition of AB field and an additional field 
described below. Potentials ( |1.6|) of the latter field are given as 
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,i;, u = D(h ,44 1} 



c^ 



ch 



r m 



A (r) 



COS</2 , 



(5.3) 



where A(r) , B (r) are arbitrary functions of r. This field is an analog of the field ( |4.1| ) in 
3 + 1 dimensions. Potentials A^\ A\ , A^ of AB field in 2 + 1 dimensions are still given 



by the formulas fll.l|). The operators 



ar dip 2, 



(5.4) 



are integrals of motions in the external field under consideration. Thus, we can impose the 
following conditions on the spinor \1/ 

(5.5) 



Po q = hk q, j 3 * = n[i-i - - ] * 



A solution of the equations Q5.1] ), ( |5.5|) has the form 



-ir 



V 



/ 



e - *^! (r) 
#2 (r) 

where the functions ipk (r) (A; = 1, 2) obey the equations 

I - 1 + A(r 



, r = k x° -(i-i Q )(p, 



-■01 (r) - (fc + 5 (r) + m) -02 (0 



, I + A (r) 
i>2 ( r ) = {ko + B (r) - m) ipi (r) 2 (r) 



(5.6) 



(5.7) 



Explicit solutions of these equations can be found in three particular cases: 

1. A (r) = B (r) = ; 

2 

2. A(r) =P = 2 ^, 7>0, fl(r) = 0; 

3. A (r) = 7T, B(r) = -, 7 > ; 

r 



(5.8) 
(5.9) 
(5.10) 



Below we consider each case in detail. 

The case ( |5.8|) corresponds to the pure AB field. For I 7^ 0, there exist bounded solutions 
of the form 



^1 (r) = ^Jk + m J| g _i| (kr) , k = ^Jk 2 -m 2 , 



1^2 (r) = eyk + m J\ s \(kr) , e = sign/ 



(5.11) 
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Here Ji(x) are Bessel functions. Solutions ( 5.11|) vanish at r = 0. They are orthogonal and 
normalized. 

For I — 0, p ^ 0, bounded solutions do not exist. However, there are unbounded at r = 
solutions. Some of them have the form 

V>i if) = \fk + m [ci J M _i (Ax) + c 2 Ji-,, (Ax)] , 

V> (r) = ^k ~m[ Cl J^ (Ax) - c 2 J^ (Ax)] . (5.12) 

Here c\, c 2 are arbitrary constants. The solutions ( |5.12j ) are still orthogonal and normalized. 
Another set of unbounded solutions (they are expressed via the Macdonald functions) reads 



^i ( r ) = ym + koKi-^ [Jm 2 — kfir ) , Re Jm 2 — k 2 , > 



ifa(r) =-y l m- k K^ y^Jm 2 - k 2 rj . (5.13) 

As we see they are defined even for some complex ko- It is interesting to remark that the 
latter solutions have finite norms. 

The case ( j5.9|) corresponds a combination of AB and uniform constant magnetic fields. 
Here we can introduce operators a£, a\~ (k = 1, 2) by relations (|2.33| ). Using the substitution 

^ ( x ) = e - ik ° x ° (p, <p) , ^2 (x) = *e"^° ^ 2 (p, <p) , (5-14) 



we present fl5.2|) in the following form 

(ho - to) (p, <p) - J^ai^ (p, <p) = , (fc + m ) ^2 (p, <p) - J^iat^i (p, V 3 ) = . (5.15) 



Now we can use the functions ( |2.47|) and the relations (]2.41|) . 

Consider first states with k 2 , ^ m 2 , I ^ 0. As in 3 + 1 dimensions, these states can be 
divided in two types ( j — 1, 2), 



V>i (p, <p) = yfk + m^2i, i_i , ^2 (p, <p) = yjko-mil;^} , fcg = m 2 + 27n , (5.16) 



where n was defined in fl2.46|) . Solutions ( p.!6| ) vanish at r = 



If k 2 , 7^ m2 i I = 0, fi ^ 0, then solutions, which are formally defined by equations ( |5.16 ) 



are unbounded at r = 0. However, they are still orthogonal and normalized. For I = 0, 
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p 7^ 0, k\ = m 2 + 7 (2A + p — 1) (A is arbitrary complex), there exist another unbounded 
solutions. They have the form ( |5.6| ), with 

~~ (5-17) 



tpi (r) = (fc + m)V>A-i,i-u ^ 2 ( r ) = V 2 ^a,m (P) 



where the functions ip\ )fl (x) are defined by flS.28p , ( |6.29[ ). These solutions are orthogonal 
and normalized as well. 

Consider now states with k$ = m 2 . Suppose ko = m; then a general solutions of Eqs. 
reads 



(r) = Nfi (p) + ce^p a , ip 2 (r)=gNe 2 p a 



MP) 



e 2 p 2 



e x x ' <ix, g 



7 



2m 2 



(5.18) 



where N, c are arbitrary constants. For p, ^ the only some states with I = have a finite 
norm. Namely the states 



■01 (r) = iVp 2 e 2 J e x x ^dx, ip 2 (r) = gAT e 2 p 2 
p 



N 



T (1 + p) sin p7T 



^ 2 (r) + 0> 2 dp 



(5.19) 



\j7r[i + p^(i)-pv(i+p) + p C? 2 ] ' y 

where ip (x) is the logarithmic derivative of T— function ( [pl|, 8.360). All the above solutions 
are singular at r = 0. One can see that lim M ^o N = 0. Thus, in 2 + 1 dimensions, in pure 
magnetic field, Dirac equation does not have solutions with ko = m, in contrast to the 
corresponding 3 + 1— dimensional case. For p ^ (in the presence of AB field) such solutions 
appear even in 2 + 1 dimensions. 

Suppose k = —m. Then, the functions 



01 (r) = 0, tfj 2 (r) = [r (l - l) 



p 2 e 2 



^ (r) + i)\ dp = 1 



(5.20) 



present solutions for I < 0. One can see that ( |5.20| ) is a particular case of ( |5.16| ). For I < 0, 
the states ( |5.20| ) vanish at r = 0. For I = 0, p ^ 0, they are singular at r = 0. If p = 0, then 
this singularity disappears, however the states do not tend to zero as r — > 0. 
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Let us finally turn to the case (|5.10|) . Here it is enough to consider 6 ^ only. Indeed, 



suppose 6 = 0; then doing the change of functions (ip — > /) 

i>i ( r ) = \/ky+m f\ (r) , ip 2 (r) = y/ko-m f 2 (r) , (5.21) 



we can transform Eqs. (|5.7|) to the form (|2.78|) with &i = y A; 2 . — f^ 2 • Solutions of Eqs. 
( f2.78j ) are given by the formulas ( |2.85| )-( [2~89| ). Thus, below we consider the case 6 ^ only. 



We introduce the following notation 



9l = yj\2l-l\+2b, g-i = y/\2l -1| -26, e = ±1 = sign/. (5.22) 
Thus solutions of Eqs. ( |5.7| ) have the form 

V>i (r) = ci/ Q+SjS (y) + c 2 I a+s -i, s+i (y) , ^2 (r) = c 3 J a+SjS (y) + c 4 / Q+s _i 5 s+1 (y) , 



, 2bk - (21-1 7 

y = 2y/m? + 'y 2 -k% r, a = l + ^ lff _ 1 , 2 + 2a + = V y (5.23) 

v/m 2 + 7 2 — /c 2 , 

where 



ci = (g e - eg_ £ ) J (21 - l)k - 267 - mg x g_ x = [ 9e — e \ c 3 , 

c 2 = - (g £ + eg-e) ^(2l-l)k -2b 1 + mg 1 g_ 1 = ^ £ + ^' £ j c 4 , (5.24) 

(g £ - eg_ £ ) 



and the relations ( |6.6| )-( |6TTT|) for the Laguerre functions I n>m were used. It follows from 



( gTgg ) that both 5(1 and are real and positive for (21 - l) > 46 2 . If (2/ - l) < 46 2 



then one of these quantities is real and positive and another one is imaginary. Suppose 
21 — l) > 46 2 . Then for the energies fc 2 < m 2 + 7 2 there exist bound states. In the latter 



case s is positive integer and k% is quantized according to the last equation ( |5.23|) . 
Another observation: For any b there exist a number I such that 

1 > g \ g \ = (21- l) 2 - 46 2 > . (5.25) 

For example, 6 = corresponds to I = 0. Then one can chose two values for a, 

a 1 = l + gig- 1 , a 2 = l-gxg-i. (5.26) 



In both cases the solutions (|5.23|) have singularity at r = but still have finite norms. 
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VI. APPENDIX 



l.The Laguerre functions J n m (x) are defined by the relation 



r(l + n) exp(-x/2) «=m,, , , N 

x 2 m, n — m + 1; x). 



(6.1) 



\\ T(l + m) r(l + n-m)' 
Here $ (a, 7; x) is the confluent hypergeometric function in a standard definition (see 
9.210). For 7 7^ —s, where s is integer and non-negative, the latter function can be presented 
by a series 



(6.2) 



to (7)* *l r(a) fe f T(7 + k) k\ 
This series is converges for any complex x. For any complex a, the Pochhammer symbols 
(ce)k are defined as follows 

T{a + k) 



(a)k = Oi{a + l)...(a + fc — 1) 



Via] 



(6.3) 



2. Let m be a non-negative integer number; then the Laguerre functions are related to 
Laguerre polynomials L%(x) ( |18|], 8.970, 8.972.1) by the equation 



In.m (^0 



r(i + m) 



\ r(i + 



n 



exp(-x/2)x^L n m - m (x), m = 0,1,2,... 



(6.4) 



1 Jn n 

g^^-a.z e~ x x n+a = 

"•' fc=0 



n! 



/ \ 

n + a 


{-x) k 


( \ 

n + a 


— k j 


k\ 


K n ) 



$(-n, 1 + a; 2) . (6.5) 



Here n are non-negative integer numbers such that 

(A 



a 



r l + a 



a(a — !)...(« — n + 1) 



r(l + n)r(l + a-n) 



n! 
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3. Using well-known properties of the confluent hypergeometric function ( ||18|| , 9.212; 
9.213; 9.216), one can easily get both relations for the Laguerre functions 



2y/x(n + l)I n+1:Tn (x) = (n-m + x)L n , m (x) - 2xl' nm {x) , (6.6) 



2^Jx(m + l)/ n , m+ i(a;) = (n - m - x)I n ^ m (x) + 2xl' n m (x) , (6.7) 

2y/xnl n ^ ljm (x) = (n - m + x)I n , m {x) + 2xV n m {x) , (6.8) 

2^/xmI n)m -i(x) = (n - m - x)I njm (x) - 2xl' n m (x) , (6.9) 

2v^m/„-i, m -i(x) = (n + m - x)I njm (x) - 2xl' n m (x) , (6.10) 
2yJ (n + l)(m + 1) I n+ljm+1 (x) = (n + m + 2 - x)I n>m (x) 

+2xl' ntm (x) , (6.11) 

and a differential equation for these functions 

Ax 2 I'l m (x) + AxI'^Jx) - [x 2 - 2x{l + n + m) + (n- m) 2 ]/ n , m (x) = . (6.12) 

Suppose I n>m (x) and I m>n (x) are linearly independent. Then, a general solution I of this 
equation has the form / = AI n ^ m (x) + BI m ^ n (x). However, whenever the condition (|6.18 



holds, I n ,m{x) and I m ,n{x) are dependent. The formulas ( |6.6j )- (|6TTT| ) and the equation (|6.12|) 



are valid for any complex n,m,x . One has to be careful applying the formulas ( |6.8|) - (|6.10| ) 
for n, m — . A straightforward calculation, which uses ( |6.1| ) and ( |6.2| ), gives 

, — , . Sill 7TITX j ~ 7" 1-fm , 

lim V n 4~im(3 ; ) = \ T(l + m)x 2 exp(x/2), lim \/mI n m _i(x) = . (6.13) 

n— >0 ' 7T v m— >0 



A combination of Eqs. ( |6 . 6| ) - ( |5T9j) results in the following relations 



2^fxl' nm (x) = \/nin-i,m (x) - Vn + l/ n +l,m( 



X 



,m— 1 (*^ ) , (6-14) 



\Jx(n+ l)/„ + i, m (x) - (n - m + x) J n>m (x) + y/xnI n -i lTn (x) = , (6.15) 
yx (m + l)/ n , m+ i (x) - (n - m - x) J niTO (x) + (x) = . (6.16) 

4. Using properties of the confluent hypergeometric function, one can get a representation 
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T(l+n) exp(x/2) 

\ tVi -l rn 1 $ 1 + n, 1 + n - m; -x , 6.17 

\ T(l + m) T(l + n — m) 



and a relation ( |]|, 9.214) 

In, m {x) = (-l) n ~ m I m , n (x), n ~ m integer . (6.18) 

5. An asymptotic formula takes place 

$(a,c;x) « ^rle x x a ~ c , Rex -> oo . (6.19) 
r(a) 

Thus we obtain the following asymptotic behavior of I n ,m{x) whenever m is not integer 

. . sinm7r /— r n+m+2 , . . „ . 

I nm (x) = Jr(l + n)r(l + m)x ~ exp(x/2), Rex^oo, (6.20) 

7T V 



and 



/„,„(.'•) =(-lP f 2 6XP{ X,2) , Rex^oo, (6.21) 
T(l + n)r(l + m) 



whenever m is integer. 

6. One can prove the following asymptotic formula 

where J v [x] are Bessel functions. 

7. Taking into account ( |6.20|) and (|6.21| ), one can see that only the functions I a+n>n (x 



lim/p +aiP+/3 I — I = Ja-p(x) , (6.22) 



with non-negative integer n and a > — 1 are quadratically integrable on the interval (0, oo) 
Such functions obey the orthonormality relation 



DC 



J Ia+n,rii^)Ia+m,m{x)dx — $m,n j (6.23) 




which follows from the corresponding properties of the Laguerre polynomials ( JT8[, 7.414.3). 
In such a case, the relation 

I a+nin (x) = ]l T ^ n + a + jj e 2 x^(x) (6.24) 
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follows from ( |6.4|) 

8. Consider a class of functions, which are closely related to Laguerre functions, and 
which appear often in various problems of mathematical physics. 

As it follows from Eq. ( |6.12| ), the Laguerre functions are solutions of the following 
eigenvalue problem 

t)c d d , . 

R a ip = Xip, R a = - — I : x— -, < x < oo, a = const. (6.25) 

Ax 4 ax ax 2 

A general solution of this problem has the form 

i/} (x) = aln tm (x) + bl m)7l (x) , a = n — m, 2A = n + m + 1, (6.26) 

where a, b are arbitrary constants. In the general case functions ip i x ) vanish as x — > oo 
only if one of the numbers n or m is positive and integer. However, one can provide such a 
behavior for any n, m, choosing some special values of a, b. Consider the functions 

1p\,a{x) = X'^W X ^{x) , 1p X ,a 0) = V>A -a (x) , (6.27) 

where W\^{x) are Whittaker functions ( JIB] , 9.220.4). The functions if)\ j(X (x) can be 
expressed via the confluent hypergeometric functions 

T (—a) xf „ /l + a 



ip\, a (x) = e 2 



I (—a) X2 ( 1 + a , \ 

—r X,l + a:x) 

r -A) V 2 ' J 



r (a) x 2 (\ — a 



r(i±2-A) V 2 

or, using ( |6.1| ), via the Laguerre functions 



— A, 1 — a; x^j 



(6.28) 



T(l + n)r(l + m) 

i>X,a \x) = — r (smn7r/ n>m (x) - smmirlr^n (x)) , 

sm (n — m)7i 

1 — Ol 1 + 01 
a = n — m, 2\ = l + n + m, n = X — , m = X — . (6.29) 

By the help of (|6.6| )-(4.16), the following properties of the functions ip\,a (x) can be estab- 
lished, 
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lp\,a 0) = V^A-i.a-l (x) + — ^ — V>A-l,a 0*0 , (6-30) 

ip x ,a (x) = v^^A-i.a+i 0*0 + - — " 2 2 V'A-l.a (a) , (6-31) 

2a^A,a (*) = (2A - 1 - x) Va !Q (a?) + - (2A - 1 - a) (2A - 1 + a) ^ A -i,a (a) , (6.32) 

2xt/> AjQ (x) = (a - x) ip x , a (x) + (2A - 1 - a) V^a-I, a+i ( x ) ( 6 - 33 ) 

= (x - 2A - 1) ^A, a - 2^A+i,a • (6-34) 

As a consequence of these properties we get 

A a 1p X ,a (X) = 2A g + " ^A-l^-l 0*0 > A^A-^a-l 0*0 = ^A,a (x) , 

A Q = — — + v^— , = ——= v^— . 6.35 

2a/x ax 2-^/x ax 

The operator R a can be expressed via the operators A a , A„, 

1 — Cfe 

= A^A a + R a ^ = A a Al - - . (6.36) 



Since ( |6.29|) is a particular case of (|6.26|) , then ij)x,a ( x ) are a l so an eigenfunctions for the 



operator R a . 

Using well-known asymptotics of the Whittaker function ( [Tj|, 9.227), we get 

\i _x r (lorl) _M 
ip\,a{x)~x 2 e a ( x -> oo; ^a,« (x) ~ - ^-x 2, a t^O, x~0. (6.37) 



2 

The functions ?/>a,o (#) have a logarithmic singularity at x ~ 0. It is important to stress 
that the functions ip\,a (x) are correctly defined and infinitely differentiable for < x < oo 
and for any complex A, a. In this respect one can mention that the Laguerre functions are 
not defined for negative integer n, m. In particular cases, when one of the numbers n, m 
is non-negative and integer, the functions ipx,a (x) coincides (up to a constant factor) with 
Laguerre functions. Thus, ip\,a {x) are eigenfunctions (of the operator R a ), which vanish at 
x — > oo . 

According to ( |6.37| ), the functions ipx,a (x) are quadratically integrable on the interval 
< x < oo whenever |a| < 1. It is not true for |a| > 1. The corresponding integrals can be 



calculated ( fll8f , 7.611), 
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r 



1 - a- 2A /X 



(A' — A) sin air 



r 



\ip x>a (x) | 2 dx 
|-0a,o (») I 2 ^ : 



r 



7T 



1 + a - 2A 



2 

l+a-2A 



/ l + a-2A' 

I CK I < 1, 



r 



a - 2A 



l-a-2A 



sin an p ( 1+Q 2 ~ 2A ) T ^-"-^ j 
^ (I - A ) 



\a\ < 1, 



X 



dx = T (1 + n) 



(6.38) 
(6.39) 
(6.40) 



Here ip(x) is the logarithmic derivative of the T— function ( ||18|| , 8.360). 

For | a | > 1, the situation is the following: the only quadratically integrable eigenfunc- 
tions of the operator R a are Laguerre functions, they also form a complete set. The functions 
^A,a {x) are orthogonal whenever arguments of the T— function in ( |6.38| ) are integer and neg- 
ative. That corresponds to n, m integer and non-negative. Thus, that is again the case 
of Laguerre functions according to (|6.29|) . If \a\ < 1, then, in the general case, the func- 
tions if)x, a (x) and ip\i, a (x), A' ^ A, are not orthogonal, as it follows from ( |6.38|) . That is a 
reflection of the fact that R a is not anymore self-conjugate operator for such values of a. 
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